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A Possible Cepheid-Like Luminosity Estimator for the Long Gamma-Ray 

Bursts 

Daniel E. Reichart^'^'^, Donald Q. Lamb^, Edward E. Fenimore^, Enrico Ramirez-Ruiz^, Thomas 

L. Cline^, and Kevin Hurley'' 

ABSTRACT 

We present a possible Cepheid-like luminosity estimator for the long gamma-ray 
bursts based on the variability of their light curves. To construct the luminosity 
estimator, we use CGi^O/BATSE data for 13 bursts, Wind/KONVS data for 5 bursts, 
Ulysses/ GRB data for 1 burst, and NEAR/XGRS data for 1 burst. Spectroscopic 
redshifts, peak fluxes, and high resolution light curves are available for 11 of these 
bursts; partial information is available for the remaining 9 bursts. We find that the 
isotropic-equivalent peak luminosities L of these bursts positively correlate with a 
rigorously-constructed measure V of the variability of their light curves. We fit a 
model to these data that accommodates both intrinsic scatter (statistical variance) and 
extrinsic scatter (sample variance). We find that L ^ V ' -o-^. If one excludes GRB 
980425 from the fit on the grounds that its association with SN 1998bw at a redshift 
of z = 0.0085 is not secure, the luminosity estimator spans ~ 2.5 orders of magnitude 
in L, and the slope of the correlation between L and V is positive with a probability of 
1 - 1.4 X 10-"^ (3.8 a). Although GRB 980425 is excluded from this fit, its L and V 
values are consistent with the fitted model, which suggests that GRB 980425 may well 
be associated with SN 1998bw, and that GRB 980425 and the cosmological bursts may 
share a common physical origin. If one includes GRB 980425 in the fit, the luminosity 
estimator spans ~ 6.3 orders of magnitude in L, and the slope of the correlation is 
positive with a probability of 1 — 9.3 x 10~^ (4.9 a). In either case, the luminosity 
estimator yields best-estimate luminosities that are accurate to a factor of ~ 4, or 
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best-estimate luminosity distances that are accurate to a factor of ^ 2. Independently 

of whether or not GRB 980425 should be included in the fit, its light curve is unique 
in that it is much less variable than the other ~ 17 light curves of bursts in our sample 
for which the signal-to-noise is reasonably good. 

Subject headings: gamma-rays: bursts 

1. Introduction 

Since gamma-ray bursts (GRBs) were first discovered (Klebesadel, Strong, & Olson 1973), 
thousands of bursts have been detected by a wide variety of instruments, most notably, the 
Burst and Transient Source Experiment (BATSE) on the Compton Gamma-Ray Observatory 
(CGRO), which will have detected about 2700 bursts by the end of CGRO^s more than 9 year 
mission in 2000 June (see, e.g., Paciesas et al. 1999). However, the distance scale of the bursts 
remained uncertain until 1997, when BeppoSAX began localizing long bursts to a few arcminutes 
on the sky, and distributing the locations to observers within hours of the bursts. This led to 
the discovery of X-ray (Costa et al. 1997), optical (van Paradijs et al. 1997), and radio (Frail et 
al. 1997) afterglows, as well as host galaxies (Sahu et al. 1997). Subsequent observations led to 
the spectroscopic determination of burst redshifts, using absorption lines in the spectra of the 
afterglows (see, e.g., Metzger et al. 1997), and emission lines in the spectra of the host galaxies 
(see, e.g., Kulkarni et al. 1998a). To date, redshifts have been measured for 13 bursts. 

Recently, Stern, Poutanen, &; Svensson (1999; see also Stern, Svensson, Sz Poutanen 1997), 
Norris, Marani 8z Bonnell (2000; see also Norris, Marani 8z Bonnell 1999), and Fenimore & 
Ramirez-Ruiz (2000; see also Ramirez-Ruiz & Fenimore 1999) have proposed trends between 
burst luminosity and quantities that can be measured directly from burst light curves, for the 
long bursts. Using 1310 BATSE bursts for which peak fluxes and high resolution light curves 
were available, Stern, Poutanen, &; Svensson (1999) have suggested that simple bursts (bursts 
dominated by a single, smooth pulse) are less luminous than complex bursts (bursts consisting of 
overlapping pulses); however, see §5. Using a sample of 7 BATSE bursts for which spectroscopic 
redshifts, peak fluxes, and high resolution light curves were available, Norris, Marani & Bonnell 
(2000) have suggested that more luminous bursts have shorter spectral lags (the interval of time 
between the peak of the light curve in different energy bands). Using the same 7 bursts, Fenimore 
Sz Ramirez-Ruiz (2000) have suggested that more luminous bursts have more variable light curves. 
These trends between luminosity and quantities that can be measured directly from light curves 
raise the exciting possibility that luminosities, and hence luminosity distances, might be inferred 
for the long bursts from their light curves alone. 

In this paper, we present a possible luminosity estimator for the long bursts, the construction 
of which was motivated by the work of Ramirez-Ruiz &; Fenimore (1999) and Fenimore Sz 
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Ramirez-Ruiz (2000). We term the luminosity estimator "Cepheid-like" in that it can be used 
to infer luminosities and luminosity distances for the long bursts from the variabilities of their 
light curves alone. We apply this luminosity estimator to every long burst in the current BATSE 
catalog in an upcoming paper (Reichart et al. 2000). 

Wc rigorously construct our measure V of the variability of a burst light curve in §2. 
Qualitatively, V is computed by taking the difference of the light curve and a smoothed version 
of the light curve, squaring this difference, summing the squared difference over time intervals, 
and appropriately normalizing the result. Our variability measure differs from that of Fenimore 
Sz Ramirez-Ruiz (2000) in three important ways: (1) we define the timescale over which the light 
curve is smoothed differently; (2) we subtract out the contribution to the variability due to Poisson 
noise differently; and (3) we combine variability measurements of light curves acquired in different 
energy bands into a single measurement of a burst's variability differently. 

We find that only smoothing timcscalcs that are proportional to burst duration appear to 
lead to significant correlations between the isotropic-equivalent peak photon luminosity L of a 
burst and V; in particular, smoothing timescales that are a fixed duration in the source frame do 
not. We take the smoothing timescale of a burst light curve, acquired in observer-frame energy 
band E, to be the time spanned by the brightest 100/% of the total counts above background, 
where / is a number between and 1 that we fix using the data in §4. We schematically illustrate 
this measure of duration in Figures 1 and 2. In Figure 1, the 50% smoothing timescale of the light 
curve is given by T^Lq 5 = Ti + T2, and the 90% smoothing timescale of the light curve is given by 
TjLq g = T3; in Figure 2, the 50% smoothing timescale of the light curve is given by TjLq = T4, 
and the 90% smoothing timescale of the light curve is given by TjLq g = T5 + Tq. We have chosen 
this measure of the smoothing timescale over others because it is robust: small variations in either 
the light curve or the value of / result in only small variations in the value of the smoothing 
timescale. This is not the case with measures like T^o and Tgo (see, e.g., Paciesas et al. 1999). 

For example, consider the case of a burst with a precursor. The value of, say, T90 can differ 
considerable if the precursor's counts above background is < 5% versus > 5% of the total counts 
above background. Likewise, if the precursor's counts above background is, say, 5% of the total 
counts above background, the duration that one measures can differ considerably if one uses Too 
versus r>9o. The effect of using a smoothing timescale that is artificially long (short) is measuring 
a variability that is artificially high (low). Our measure of the smoothing timescale does not suffer 
from such robustness problems. 

We present our measure of the isotropic-equivalent peak photon luminosity L of a burst in 

§3. In §4, we expand the original Ramirez-Ruiz & Fenimore (1999) sample of 7 bursts to include 
a total of 20 bursts, including 13 BATSE bursts, 5 Wind/KONVS bursts, 1 Ulysses/GRB burst, 
and 1 NEAR/XGRS burst. Spectroscopic redshifts, peak fluxes, and high resolution light curves 
are available for 11 of these bursts; partial information is available for the remaining 9 bursts. 

Also in §4, we construct our luminosity estimator, which differs from that of Fenimore &i 
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Ramirez-Ruiz (2000) in two important ways: (1) applying the Bayesian inference formalism 
developed by Reichart (2000), we fit a model to the data that accommodates both intrinsic scatter 
(statistical variance) in two dimensions and extrinsic scatter (sample variance) in two dimensions; 
and (2) again applying this Bayesian inference formalism, we determine not only the best estimate 
for L as a function of V, but also the uncertainty in L as a function of V, as well as approximate 
these functions with analytic expressions. We state our conclusions in §5. 



2. The Variability Measure 

We now rigorously construct a measure of the variability of a burst light curve. We require it 
to have the following properties: (1) we define it in terms of physical, source-frame quantities, as 
opposed to measured, observer-frame quantities; (2) when converted to observer-frame quantities, 
all strong dependences on redshift and other difficult or impossible to measure quantities cancel 
out; (3) it is not biased by instrumental binning of the light curve, despite cosmological time 
dilation and the narrowing of the light curve's temporal substructure at higher energies (Fenimore 
et al. 1995); (4) it is not biased by Poisson noise, and consequently can be applied to faint bursts; 
and (5) it is robust; i.e., similar light curves always yield similar variabilities. We also derive 
an expression for the statistical uncertainty in a light curve's measured variability. Finally, we 
describe how we combine variability measurements of light curves acquired in different energy 
bands into a single measurement of a burst's variability. 

We first define the variability of a burst light curve, acquired in observer-frame energy band 
E, in terms of physical, source-frame quantities: 



where L^^^tg) is the luminosity of the burst in source-frame energy band Eg = E{1 + z) as a 
function of source-frame time tg, and (L^' ikSf){tg) is the convolution of this function and a boxcar 
smoothing function of area equal to 1 and width equal to the smoothing timescale of L^''{tg), i.e., 
the source-frame smoothing timescale of the light curve. Since our definition of the smoothing 
timescale is a robust measure of duration (§1), Equation ([l|) is a robust measure of variability. 

Next, we convert Equation ([l|) to observer-frame quantities: 

^'"^ jZo[c''{to)-B^{to)nAnD^{z){i+zr]^{i+z)^dto ■ 

Here, C^{to) and B^{to) are the total (source plus background) and background photon count 
rates in observer-frame energy band ii^ as a function of observer- frame time to, [{C^ — B^)'k Sf]{to) 
is the convolution of C^{to) — B^(to) and a boxcar smoothing function of area equal to 1 and 
width equal to the smoothing timescale of C^{to) — B^{to), i.e., T^, AQ is the effective solid 
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angle of the emitted light, D(z) is the comoving distance to the burst at redshift z, a is either 1 or 
2 depending on whether L^''{ts) is a photon number or energy luminosity, and /? ~ 1 — 0.4 = 0.6: 
the factor of (1 + zy is due to cosmological time dilation, and the factor of « (1 + z)~^-'^ is due 
to the narrowing of the light curve's temporal substructure at higher energies (Fenimore et al. 
1995). By construction, the strong dependences on Af2, z, a, and /?, and a weak dependence on 
the cosmological model, cancel out. Furthermore, since B^{to) always varies over much longer 
timescales than Tj', Equation (|2|) simplifies to: 

..E _ rooiC' {to) -{C^^Sf){to)f dtp 

We distinguish Equations @) and dH) from Equation (|^) with a subscript P because they differ 
in one fundamental way: unlike L^^{ts) in Equation (|l|), C^{to) suffers from Poisson noise. This 
additional source of variability biases Equations (||) and @ to higher values. Although this bias is 
minor for bursts that are significantly brighter than the background, it cannot be ignored for faint 
bursts. We determine and subtract out the contribution to the variability due to Poisson noise 
below. 

However, before we can subtract out the contribution to the variability due to Poisson noise, 
we must first bin the light curve, since observed light curves typically are discrete, not continuous. 
Consequently, we replace the integrals in Equation (^) with sums from i = lto N, C^iQ with 
Cf , B^{to) with Bf, and {C^*Sf){to) with its discrete equivalent, Si{Cf,Nf). Here, Cf is the 
light curve to be smoothed, N^" is the smoothing timescale, i.e., T^, but measured in number of 
bins, and 



5.(Cf,iV.) = ^ 



j=i-n 



(4) 



where Ux is the truncated integer value of {Nx — l)/2. Binning can bias the variability measure, 
since binning wipes out all variability on timescales shorter than the sampling timescale of the 
light curve, and the effective sampling timescale of the light curve in the source frame decreases 
with increasing burst redshift: although the sampling timescale of the light curve is fixed by 
the instrument in the observer frame, (1) cosmological time dilation and (2) the narrowing of 
the light curve's temporal substructure at higher energies (Fenimore et al. 1995) decrease the 
effective sampling timescale of the light curve in the source frame. We remove this redshift bias by 
smoothing all light curves on a source-frame timescale of 64 ms (the shortest sampling timescale of 
BATSE), which corresponds to a timescale of 64(1 + z)^ ms in the observer frame. Consequently, 
we replace Cf' with Si{Cj', Nz), where Nz = 64(1 + z)^ /At, and At is the observer-frame sampling 
timescale of the light curve in milliseconds. Equation (H) becomes: 



^ J:ll[S^iCf,Nz)-S.iCf,Nf )]^ 
^'"^ Ell[S^{Cf,Nz)-Bf]^ 
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where we take Cf and to be measured in counts, not counts per unit time. Although the 



computation of Vpp now depends on z and P, these dependences are very weak: wide variations in 



the values of these parameters do not significantly change the measured variabilities of the bursts 
in our sample; i.e., burst light curves have very little power on such short timescales. Likewise, 
wide variations in the value of the effective sampling timescale of the light curve that we impose 
in the source frame also have little effect on the measured variabilities of these bursts. 

Now, we determine and subtract out the contribution of Poisson noise to the variability, and 
simultaneously we derive an expression for the statistical uncertainty in a light curve's measured 
variability. First, we rewrite the expressions Si{Cj' , Nz) — Si{Cj",Nj") and Si{Cj' , Nz) — Bf from 
Equation ^ as weighted sums of the statistically independent measurements Cj: 

'^"^ ~ ElAEUh^c, - B^ ^ ^ 

where aij and bij are weights that differ for each burst, but can be computed straightforwardly using 
Equation (^). Since we are in the Gaussian limit, the uncertainty in each Cj is simply it-y/Cj. Since 
each Cj is statistically independent, their weighted uncertainties can be summed in quadrature. 

N „ ri ^ /^N 



Hence, the sums over j in Equation (|6|) and their uncertainties are: J2j=i o-ijCj ib yJ2j=i '^ijCj 



and Ef=i h,C, -B,± y'Ef=i h%C,. 

Next, in accordance with Equation (|6|), we square these expressions, which yields expressions 
consisting of three terms: the square of the original term, a positive term due to Poisson noise, 
and an uncertainty term. Subtracting out the contributions to Equation (^) due to Poisson noise 
yields: 



_ E^Ii[(I]f=i o-ijCj)'^ - J2jLi (ifjCj] 
' " Ef=i [(Ef=i b.,c, - B,y - Ef=i yi,c,] ' ^ ^ 



where the uncertainty in the ith term of the numerator is ±2{J2j=i 0'ijCj)yJ2j=i ^ij^ji ^'^^ the 

uncertainty in the ith term of the denominator is ±2(X]^i i>ijCj — Bi)^J2jLi ^fj Cj • Recognizing 
that only ~ N/Nz of the terms in the sums over i are statistically independent, since the sums 
over j correspond to a convolution of the light curve and a boxcar smoothing function of width 
equal to Nz bins, the uncertainties in the sums over i are a factor of \fNz larger than what 
would be derived if all of the terms were statistically independent. Hence, the uncertainty in the 
numerator is ±2^^ Nz Ei^i(EjLi ^.ijCjY Yl!j=i o-ljCj, and the uncertainty in the denominator is 

±2^ Nz J2iLi{J2jLi ^ijCj — BiY EjLi bfjCj. These uncertainties can be at most only weakly 
correlated. Taking them to be independent, we find that the statistical uncertainty in a light 
curve's measured variability is approximately given by: 



(TyE = VF 



4iV. EiIi(Ef=i a^jC.r Ef=i a2,C, ^ 4iV, EiIi(Ef=i h.c, - B,Y Y.U 



\ {EiIi[(Ef=i a.,c,Y - Ef=i alc,]Y {EiIi[(Ef=i h,c, - b,y - Ef=i hlc,]Y ' 

(8) 
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All that remains is to describe how we combine variability measurements of Hght curves 
acquired in different energy bands, V^, into a single measurement of a burst's variability: Vj. 
For each burst in our sample, light curves were acquired in typically 3 or 4 independent energy 
bands (see Table 2). We find that (1) the smoothing timescales of these bursts' light curves 
decrease with energy as ~ E~^'^ (Figure 3), and (2) the variabilities of these bursts' light curves 
are approximately constant across energy bands (Figure 4). In hindsight, the former result is not 
surprising, given our definition of the smoothing timescale and the principle result of Fenimore 
et al. (1995); however, this does constitute an independent confirmation of their result. The 
latter result suggests that we can model a burst's variability as a constant across energy bands. 
Applying the Bayesian inference formalism developed by Reichart (2000) for fitting models with 
extrinsic scatter (sample variance) to data with intrinsic scatter (statistical variance, in this case 
given by Equation |8|), we fit this model to the typically 3 or 4 independent measurements of a 
burst's variability made in independent energy bands, Vj' , resulting in a single measurement of 
that burst's variability, Vf, and the uncertainty in Vf. We plot the > 25 keV light curves of the 
most and least variable cosmological BATSE bursts in our sample in Figure 5. 



3. The Luminosity Measure 

Let P be the peak flux of a burst in photons cm~^ s~^ between observer-frame energies Ei and 
Eu. The isotropic-equivalent peak photon luminosity of the burst in erg between source-frame 
energies 100 and 1000 keV is given by 



J^;HE)dE 



(9) 



where $(£') is the observer-frame spectral shape, which we parameterize with the Band function 
(Band et al. 1993). For each burst in our sample, we compute the value of L for each of the 
54 parameterizations of the Band function in Band et al. (1993); the luminosity error bars in 
Figures 9 and 10 are dominated by this variation in the parameterization of the Band function, 
but clearly are too small to matter. Hence, L is very insensitive to reasonable variations in the 
parameterization of the Band function for the bursts in our sample. The choice of source-frame 
energy range also is very unimportant: we chose 100 - 1000 keV to approximately match the 
observer-frame energy range in which BATSE measures peak fluxes, 50 - 300 keV, for redshifts 
that are typical of the bursts in our sample: z ~ 1 — 2. 



4. The Luminosity Estimator 



We list our sample of 20 bursts in Table 1; it consists of every burst for which redshift 
information is currently available. Peak fluxes are available for all 20 bursts. Spectroscopic 
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redshifts and 64-ms or better resolution light curves are available for 11 bursts. Spectroscopic 
redshifts but only 1-s resolution light curves are available for 2 bursts. We compute variability 
lower limits for these bursts: we compute variabilities from their 1-s resolution light curves without 
further degrading the effective resolution of these light curves by smoothing them on the 64-ms 
source- frame timescale (§2); these variabilities are lower limits to the variabilities that we would 
compute if 64-ms or better resolution light curves were available. Redshift upper limits (1) from 
the non-detection of the Lya forest in host galaxy spectroscopy in the case of 1 burst, and (2) 
from the non-detection of the Lyman limit in afterglow and host galaxy photometry in the case 
of 6 bursts, and 64-ms or better resolution light curves are available for the remaining 7 bursts. 
We compute luminosity upper limits for these bursts. We compute variabilities for these bursts 
for both z = and z equal to the redshift upper limit; for all 7 bursts, both of these variabilities 
are nearly identical, testifying to the weakness of the computational dependence of the variability 
on redshift (§2). In addition to the data listed in Table 1, we used Ulysses/GHB data for GRBs 
970228, 990712, and 991208 to test the consistency of our results across instruments; we find the 
measured variabilities and luminosities of these bursts to be fully consistent across instruments. 

We now construct the luminosity estimator. First, we compute the variabilities, Vj, and 
the isotropic-equivalent peak luminosities, L,0 of the bursts in our sample. We compute the 
variabilities as functions of /, where we use the value of / to compute the smoothing timescales, 
, of all of the light curves in our sample (§1). We identify values of / that lead to successful 
and robust luminosity estimators below. Regardless of the value of / that we use to compute the 
variabilities, the distribution of our sample's bursts in the logL-logVy plane appears to be well 
modeled by a normal distribution about a straight line (see, e.g.. Figure 9). We parameterize the 
line by: 

log Vf{L) =logVf + b + m{logL -logL), (10) 

where b is the intercept of the line, m is its slope, and Vf and L are the median values of Vj and 
L for the bursts in our sample for which spectroscopic redshifts, peak fluxes, and 64-ms or better 
resolution light curves are available. We parameterize this line as a function of L, instead of as a 
function of Vf, because the data do not fully rule out the possibility of m ^ 0, whereas the data do 
rule out the possibility of very large values of m. We parameterize the normal distribution about 
this line by aiog Vf j which is half of the distribution's l-cr width along the log Vf axis. Applying the 
Bayesian inference formalism developed by Reichart (2000) for fitting data with extrinsic scatter 
(in this case, CT\ogVf) to data with intrinsic scatter (in this case, the uncertainties in the measured 
values of Vf and L), we determine values and uncertainties for the model parameters (6, m, and 
o'logvy) as functions of / (Figure 6).^ 



*We use peak fluxes measured on a 1-s timescale, and we take Ho — 65 km s ^ Mpc ^, 57m ~ 0.3, and Qa = 0.7; 
the luminosity estimator is very insensitive to these choices. 

^This Bayesian inference formalism deals only with measurements with Gaussian error distributions, not with 
lower or upper limits. However, this formalism can be straightforwardly generalized to deal with limits as well, using 
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Of the 20 bursts in our sample, GRB 980425 is unique due to its possible association with 
SN 1998bw (Kulkarni et al. 1998b; Galama et al. 1998) at a redshift of z = 0.0085 (Tinney et al. 
1998); the 12 other bursts for which spectroscopic redshifts are available have 0.430 < z < 3.418. 
Consequently, we first construct the luminosity estimator excluding GRB 980425 from the above 
fits; we then repeat the fits including this burst. Excluding GRB 980425 from the fits, we first 
identify values of / that lead to successful and robust luminosity estimators. One measure of the 
success of a luminosity estimator is the probability that its slope, m, departs from m = 0. Since 
the best-fit values of m are positive (Figure 6), we compute and plot in Figure 7 the probability 
that m < 0, P{m < 0), as a function of /. We find that values of / between 0.2 and « 0.5 
lead to luminosity estimators with slopes that are positive with approximately equally probability: 
P{m < 0) « 1.4 X 10~^ (3.8 a). Values of / ^ 0.2 and / ^ 0.5 lead to luminosity estimators 
with slopes that are positive with increasingly less probability. In addition to having a slope that 
is positive (or negative) with large probability, a successful luminosity estimator also should be 
robust: small changes in the value of / should lead to only small changes in the estimated values 
of L. Consequently, we also plot in Figure 7 the extrinsic scatter of the data along the logL axis, 
ciogL = o"iogV//"^) as a function of /; this function is a measure of the maximum amount that 
the estimated luminosities of the bursts in our sample can change as a function of /. We find 
that values of / ^ 0.5 lead to increasingly less robust luminosity estimators. Taking / = 0.45, 
which corresponds to the minimum in P{m < 0), we compute smoothing timescales, for each 
of the instrument's available energy bands, Ei, (Table 2), and using these smoothing timescales, 
we find that that log?/ = -1.088, log! = 51.852, b = 0.013l:[j:[j^^, m = 0.302^0^^11, and 

3 3+1.1 

o"iogV7 ~ ^--'-'''^-aole- Hence, we find that L ~ Vj^^q'^^ . We plot the two-dimensional uncertainty 
distributions of the fitted values of the model parameters in Figure 8. We plot the data and the 
best-fit model of the distribution of these data in the logL-log V/ plane for / = 0.45 in Figure 9. 

Although the solid line in Figure 9 is a good approximation to the best estimate for L as 
a function of Vf, the dotted lines in Figure 9 do not correspond to the uncertainty in L as a 
function of Vf. The dotted lines in Figure 9 mark the best-fit width, given by criogv^ = 0.175, of 
the best-fit model of the distribution of the data in the log L-log Vf plane; they do not account 
for the uncertainties in the fitted values of the model parameters, and the uncertainties in the 
fitted values of b and m in particular. Applying the Bayesian inference formalism of Reichart 
(2000), we formally compute the best estimate for L as a function of Vf, and the uncertainty 
in L as a function of Vf, which we plot as the solid curves in Figure 10. The computation of 
this distribution in the log L-log Vf plane is numerically intensive, so we also provide analytic 
approximations to these functions, which we plot as the dotted curves in Figure 10. The best 
estimate for L as a function of Vf can be approximated by Equation (|lO|), using the best-fit values 
of the model parameters. Approximating the uncertainty distributions of the fitted values of the 



two facts: (1) a limit can be given by the convolution of a Gaussian distribution and a Heavyside function, and (2) 
convolution is associative. In this paper, we fit to both measurements and limits. 
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model parameters (Figure 8) as uncorrelated and Gaussian, the uncertainty in L as a function of 
Vf can be approximated by (Reichart 2000): 



(^\ogVf{L) = ^Jf^lgYf + + ^U^ogL- logL)2, (11) 

where (TiogV/(-^) is the 1-a uncertainty in Vj as a function of L, criagVf = 0.175, af, 0.084 is 
the uncertainty in the fitted value of b, and am ~ 0.094 is the uncertainty in the fitted value 
of m. At its best, the luminosity estimator constructed excluding GRB 980425 from the fits 
yields best-estimate luminosities that are accurate to a factor of ~ 4, or best-estimate luminosity 
distances that are accurate to a factor of ~ 2. 

In addition to being used to estimate the luminosities of bursts for which spectroscopic 
redshifts are not available, and the uncertainties in these luminosities. Figure 10 also can be used 
to determine if other bursts for which spectroscopic redshifts are available, such as GRB 980425 or 
future bursts, are consistent with the above-fitted model. Consequently, we also plot GRB 980425 
in Figure 10. Although its implied luminosity, L = 5.2 ± 2.0 x 10^^ erg s~^, is « 4 — 6 orders of 
magnitude less than the 12 other measured luminosities of bursts in our sample, its variability, 
Vf=QA5 = 0.00 lb 0.01, also is considerably less than the ~ 17 other reasonably well-constrained 
variabilities of bursts in our sample, and clearly is consistent with the above-fitted model. We 
plot the > 25 keV BATSE light curve of GRB 980425 in Figure 11. Repeating the fits including 
GRB 980425 yields: / = 0.47, logV) = -1.073, L = 51.697, b = -0.035l|]:^, m = 0.306l[j:[]^?, 
o-iogV7 = O-l'^O-aoIi' and P{m < 0) = 9.3 x 10"^ (4.9 a) (Figures 6 and 7). Hence, we find that 

3 3+^" 

L ~ Vj^qIy ■ We also plot the two-dimensional uncertainty distributions of the fitted values of 
the model parameters in Figure 8. We plot the data and the best-fit model of the distribution of 
these data in the logL-logV/ plane for / = 0.47 in Figure 12. We plot the best estimate for L 
as a function of Vj, the uncertainty in L as a function of Vj, and our analytic approximations to 
these functions, given by Equations ([lO|) and (11) and the fitted values of the model parameters, 
in Figure 13. 



5. Discussion and Conclusions 

We have presented a rigorously-constructed measure of the variability of a burst's light curve. 
Using this variability measure and a sample of 20 bursts, consisting of every burst for which 
redshift information is currently available, we have shown that a significant correlation exists 
between the variability of a burst's light curve, and the burst's isotropic-equivalent peak photon 
luminosity. This correlation between variability and luminosity is in agreement with the trends 
found by Stern, Poutanen, & Svensson (1999) and Fenimore & Ramirez-Ruiz (1999). That is, 
more variable ("complex" in the terminology of Stern, Poutanen, & Svensson 1999) bursts are 
more luminous, while less variable ("simple" in the terminology of Stern, Poutanen, & Svensson 
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1999) bursts are less luminous.]^ 

Furthermore, from the correlated variabilities and luminosities of our sample of 20 bursts, we 
have constructed a possible Cepheid-like luminosity estimator for the long bursts. If one excludes 
GRB 980425 from the fits, the luminosity estimator spans ~ 2.5 orders of magnitude in luminosity, 
and its slope is positive with a probability of 1 — 1.4 x 10~^ (3.8 a). GRB 980425, however, is 
consistent with the fitted model. If one includes this burst in the fits, the luminosity estimator 
spans « 6.3 orders of magnitude in luminosity, and its slope is positive with a probability of 
1 -9.3 X lO-'^ (4.9 a). 

Future bursts will either increase these probabilities, or possibly disprove the luminosity 
estimator. However, since the uncertainty in L as a function of Vf in Figures 10 and 13 is 
dominated by extrinsic scatter (i.e., criogVf)^ and not by the uncertainties in the fitted values of 
the model parameters, a larger sample of bursts is unlikely to improve the predictive power of the 
luminosity estimator. Currently, the luminosity estimator yields best-estimate luminosities that 
are accurate to a factor of ~ 4, or best-estimate luminosity distances that are accurate to a factor 
of ^ 2. 

However, independently of whether or not the luminosity estimator is eventually disproved, 
the light curve of GRB 980425 is unique in that it is much less variable than the other ~ 17 light 
curves of bursts in our sample for which the signal-to-noise is reasonably good: logyy=o.45 < —1.5 
with a probability of 1 — 3.4 X 10~^ (3.4 a), and log Vy=o.45 < —1 with a probability of 1 — 3.5 x 10^^^ 
(9.9 0"). The argument has been made that the association of GRB 980425 with SN 1998bw at 
the unusually low redshift of z = 0.0085 probably is accidental because the light curve of GRB 
980425 is no different than the light curves of the cosmological bursts. On the contrary, we find 
that the light curve of GRB 980425 is different from the light curves of the cosmological bursts. 
Consequently, GRB 980425 may well be associated with SN 1998bw. 

If GRB 980425 is associated with SN 1998bw, and if the luminosity estimator is correct, the 
fact that GRB 980425 is consistent with the fitted model suggests that GRB 980425 and the 



However, we must draw attention to a potential disagreement between the primary conclusion of Stern, Poutanen, 
& Svensson (1999), and the width of the luminosity distribution of the multiply-peaked bursts with spectroscopically- 
measured redshifts in our sample. Stern, Poutanen, & Svensson (1999) find that the differential peak count rate 
distribution of their complex, or multiply-peaked BATSE bursts peaks about a factor of 4 above threshold, while the 
differential peak count rate distribution of their simple, or singly-peaked BATSE bursts does not have a similar peak. 
They interpret this to mean that their complex bursts are more luminous and at higher redshifts than their simple 
bursts: they argue that this peak corresponds to the peak in the star-formation history of the universe at z ~ 1.5, and 
that it is not a threshold effect. However, as this peak is narrow, spanning less than an order of magnitude, this could 
only be the case if the luminosity function of the complex bursts is similarly narrow; otherwise, this feature would be 
washed out. We have visually examined the light curves of the bursts with spectroscopically-measured redshifts in 
our sample, and find that only GRB 980425 and GRB 970508 appear to be singly-peaked. The remaining, multiply- 
peaked bursts span ~ 2.5 orders of magnitude, which appears to contradict their interpretation of this peak that they 
find in the differential peak count rate distribution of their complex bursts. 
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cosmological bursts may share a common, or at least a related, physical origin, although they 
cannot share a common redshift distribution and/or luminosity function (Graziani, Lamb, & 
Marion 1999). This conclusion is made more intriguing by the recent discoveries of supernova-like 
components to the late afterglows of the cosmological bursts GRB 970228 (Reichart 1999, Galama 
et al. 2000; Reichart, Castander, k, Lamb 2000; Reichart, Lamb, k, Castander 2000) and GRB 
980326 (Bloom et al. 1999). 
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Table 1. Redshifts, Luminosities, and Variabilities 



GHB 


z 


r a 


T T 

Vf=0.45 


Instrument 


z Reierence 


970228 


0.695 


5.0 ± 1.7 X lO'^^ 


0.08 ± 0.05 


Mnd/KONUS 


Djorgovski et al. 1999b 


970508 


0.835 


1.0 ±0.2 X lO'^^ 


0.05 ± 0.02 


CGRO/BATSE 


Motzgcr ct al. 1997 


970828 


0.958 


7.1 ± 1.4 X 10^5^ 


0.10 ±0.00 


CGRO/BATSE 


Djorgovski 1999 


971214 


3.418 


5.8 ± 1.4 X 10^5^ 


0.11 ±0.01 


CGRO/BATSE 


Kulkarni et al. 1998a 


980308 


< 4.3*' 


< 7.4 X 10^^ 


0.24 ± 0.02 


CGRO/BATSE 


Schaefer et al. 1999 


980326 


< 4.3'' 


< 2.5 X lO''^ 


0.09 ±0.02 


CGRO/BATSE 


Groot et al. 1998 


980329 


< 3.9^^ 


< 5.7 X 10^^ 


0.04 ±0.01 


CGRO/BATSE 


Djorgovski 1999 


980425 


0.0085 


5.2 ± 2.0 X 10*® 


0.00 ± 0.01 


CGRO/BATSE 


Tinney et al. 1998 


980519 


< 4.3^= 


< 2.2 X 10^^ 


0.16 ±0.03 


CGRO/BATSE 


Halpcrn ot al. 1999 


980613 


1.096 


1.3 ±0.2 X 10^^ 


> 0.03^^ 


CGRO/BATSE 


Djorgovski et al. 1999a 


980703 


0.967 


3.6 ±0.7 X 10" 


0.06 ± 0.01 


CGRO/BATSE 


Djorgovski et al. 1998 


981220 


< 4.3'' 


< 4.3 X 10^^ 


0.08 ± 0.03 


Wind/KONUS 


Masctti ct al. 1998 


990123 


1.600 


8.5 ± 1.0 X 10^^ 


0.10 ±0.00 


CGRO/BATSE 


Kulkarni et al. 1999 


990510 


1.619 


4.3 ± 0.5 X 10^^ 


0.24 ±0.01 


CGRO/BATSE 


Beuermann et al. 1999 


990705 


< 5.5'^ 


< 3.5 X lO'^^ 


0.15 ±0.05 


Ulysses/GRB 


Masetti et al. 2000 


990712 


0.430 


3.8 ± 1.5 X 10^" 


0.02 ± 0.03 


Mnd/KONUS 


Galama et al. 1999 


991208 


0.706 


2.9 ± 1.0 X lO'^^ 


0.08 ± 0.03 


Wind/KONVS 


Djorgovski et al. 1999c 


991216 


1.020 


1.1 ± 0.2 X lO'^^ 


0.14 ± 0.01 


CGRO/BATSE 


Vreeswijk et al. 1999 


000131 


< 5.5*= 


< 6.7 X 10^^ 


0.11 ±0.01 


Wmd/KONUS 


Pedersen et al. 2000 


000301 


2.034 


6.2 ± 2.0 X 10^2 


> 0.03*^ 


NEAR/XGRS'' 


Castro et al. 2000 



''Isotropic-equivalent peak photon luminosity in erg between source- frame energies 100 and 
1000 keV, for peak fluxes measured on a 1-s timescale, Ho = 65 km Mpc~^, Q,m = 0.3, and 
SIa = 0.7; upper limits are 1 cr. 

''Prom non-detection of Lyman limit in afterglow or host galaxy photometry. 

"^EYom non-detection of Lya forest in host galaxy spectroscopy. 

''Only 1-s resolution light curve is available; lower limit is 1 cr. 

''Peak flux from Ulysses/GRB. 
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Table 2. Smoothing Timescales 



GRB 


Instrument^ 


rpEi b 

/=0.45 


rpE^ b 
/=0.45 


rpE^ b 
-' /=0.45 


^ /=0.45 


970228 


Mnd/KONUS 


3.235 


2.891 


1.061 


c 


970508 


CGRO/BATSE 


3.711 


3.114 


2.944 


1.333 


970828 


CGRO/BATSE 


19.057 


15.875 


12.335 


11.276 


971214 


CGRO/BATSE 


7.769 


6.706 


5.943 


2.080 


980308 


CGRO/BATSE 


4.987 


6.135 


4.903 


_d 


980326 


CGRO/BATSE 


0.782 


0.836 


_e 


_d 


980329 


CGRO/BATSE 


5.371 


5.033 


4.769 


4.382 


980425 


CGRO/BATSE 


7.819 


6.701 


5.169 


0.839 


980519 


CGRO/BATSE 


8.286 


7.232 


4.742 


2.872 


980613 


CGRO/BATSE 


8.160 


8.432 


6.528 




980703 


CGRO/BATSE 


18.024 


16.659 


15.191 


11.733 


981220 


Wind/KONVS 


2.152 


1.871 


0.740 


C 


990123 


CGRO/BATSE 


22.673 


19.346 


13.340 


7.345 


990510 


CGRO/BATSE 


6.911 


4.964 


3.118 


_d 


990705 


Ulysses/GRB 


4.986 


_f 


_f 


_f 


990712 


Mnrf/KONUS 


2.405 


3.595 


0.925 


_c 


991208 


Mnd/KONUS 


7.737 


4.682 


4.299 


_c 


991216 


CGRO/BATSE 


4.641 


3.842 


2.905 


2.009 


000131 


Wind/KONUS 


6.151 


4.632 


2.810 


C 


000301 


NEAR/XGRS 


1.552 


_f 


_f 


_f 



"CGRO/BATSE: 25 keV < Si < 55 keV, 55 keV < ^2 < 110 keV, 110 
kcV <^ Es ^ 320 keV, E4 ^ 320 keV; Mnrf/KONUS: 10 keV ^ -Ei ^ 45 
keV, 45 keV < E2 < 190 keV, 190 keV < -B3 < 770 keV; Ulysses/GRB: 
22 keV <Ei^< 15(rkeV; NEAR/XGRS: 300 keV < Ei < 1000 keV. 

''45% smoothing timescale in seconds of light curve in observer-frame 
energy band Ei. 

'^Only three energy bands are available for this instrument. 

"^The smoothing timescale of the light curve in this energy band is 
less than the imposed effective sampling timescale of the light curve due 
to negligible emission in this energy band; consequently, we drop this 
energy band. 

^There is an unusual, positive, constant, systematic offset of a portion 
of the light curve in this energy band; consequently, we drop this energy 
band. 

^Only one energy band is available for this instrument. 
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Fig. 1. — A schematic illustration of our measure of burst duration. The solid line represents a 
burst light curve, acquired in observer-frame energy band E, and the dashed line represents the 
background. Top panel: The hashed area is the brightest 50% of the total counts above background. 
Hence, the 50% duration, or smoothing timescale, of this light curve is given by T^Lq^ = Ti +T2. 
Bottom panel: The hashed area is the brightest 90% of the total counts above background. Hence, 
the 90% duration, or smoothing timescale, of this light curve is given by T^^^ g = Ta. 

Fig. 2. — Another schematic illustration of our measure of burst duration. The solid line represents 
a burst light curve, acquired in observer-frame energy band E, and the dashed line represents the 
background. Top panel: The hashed area is the brightest 50% of the total counts above background. 
Hence, the 50% duration, or smoothing timescale, of this light curve is given by7/^Q 5 = T4. Bottom 
panel: The hashed area is the brightest 90% of the total counts above background. Hence, the 90% 
duration, or smoothing timescale, of this light curve is given by T^^q 9 = T5 + Tq. 

Fig. 3. — Bursts in our sample for which spectroscopic redshifts, peak fluxes, and 64-ms or better 

resolution light curves are available. For each of these bursts, light curves were acquired in typically 
3 or 4 independent energy bands E. For each burst, we compute the durations, or smoothing 
timescales, T^Lq 45 of its 3 or 4 independent light curves, and plot 45 as a function of source- 
frame energy band Eg = E{1 + z). The horizontal bars mark the source-frame energy bands. We 
plot solid lines if the light curve has > 3000 total counts above background, and dotted lines if the 
light curve has < 3000 total counts above background. The dotted portions of these curves suffer 
from low signal-to-noise; very little variability information can be gleaned from these light curves. 
However, the trend is clear: the durations, or smoothing timescales, TjL^ 45 of these bursts are 
shorter at higher energies. From the solid portions of these curves, we find that Tp^^^^ ~ 
in agreement with the principle result of Fenimore et al. (1995). Other values of / yield similar 
results. 

Fig. 4. — Bursts in our sample for which spectroscopic redshifts, peak fluxes, and 64-ms or better 
resolution light curves are available. For each of these bursts, light curves were acquired in typically 
3 or 4 independent energy bands E. For each burst, we compute the variabilities VjL_Q 4^ of 
its 3 or 4 independent light curves, and the combined variability V/=o.45 of the burst. We plot 
log 45 — log V/=o.45 versus source-frame energy band Eg = E{\ + z). The horizontal bars mark 
the source-frame energy bands. Clearly, there is no significant trend with energy. Other values of 
/ yield similar results. 

Fig. 5.— The > 25 keV hght curves of the most (GRB 990510) and least (GRB 970508) variable 
cosmological BATSE bursts in our sample. In the case of GRB 990510 {z = 1.619), we find that 
V/=o.45 = 0.24 ± 0.01. In the case of GRB 970508 (z = 0.835), we find that V/=o.45 = 0.05 ± 0.02. 

Fig. 6. — The value of \ogVf as a function of /, and the best-fit values, and 1 cr uncertainties in 
the best-fit values, of the model parameters 6, m, and (JiogVf as functions of /. In the left panels, 
we exclude GRB 980425 from the fits; in the right panels we include GRB 980425 in the fits. 
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Fig. 7. — The probability P{m < 0) that m < as a function of /, and the extrinsic scatter (Jiogi 
of the data along the logL axis as a function of /. In the left panels, we exclude GRB 980425 from 
the fits; in the right panels we include GRB 980425 in the fits. 

Fig. 8. — The two-dimensional uncertainty distributions of the fitted values of the model parameters 
b, m, and criogv^- The curves mark the 1 a credible regions. The solid curves are for the case in 
which we exclude GRB 980425 from the fits (/ = 0.45); the dotted curves are for the case in which 
we include GRB 980425 in the fits (/ = 0.47). 

Fig. 9. — The variabilities Vf=oA5 and isotropic-equivalent peak luminosities L of the bursts in our 
sample, excluding GRB 980425. The solid and dotted lines mark the center and 1 a widths of the 
best-fit model distribution of these bursts in the log L-log y/-=o.45 plane. 

Fig. 10. — The solid lines mark the best estimate for L as a function of Vf=o.45, and the 1 a 
uncertainty in L as a function of Vy=o.45 for the case in which we exclude GRB 980425 from the 
fits. The dotted lines mark analytic approximations to these functions, given by Equations ( p!o[ ) 



and (11) and the fitted values of the model parameters. Clearly, GRB 980425 is consistent with 
the fitted model. 

Fig. 11.— The > 25 keV BATSE light curve of GRB 980425. We find that Vf=oA5 = 0.00 ± 0.01. 

Fig. 12. — The variabilities V/=o.47 and isotropic-equivalent peak luminosities L of the bursts in 
our sample, including GRB 980425. The solid and dotted lines mark the center and 1 a widths of 
the best-fit model distribution of these bursts in the log L-log V/=o.47 plane. 

Fig. 13. — The solid lines mark the best estimate for L as a function of Vj=o.47) and the 1 a 
uncertainty in L as a function of Vj=o.47 foi' the case in which we include GRB 980425 in the fits. 
The dotted lines mark analytic approximations to these functions, given by Equations (^) and 



(11) and the fitted values of the model parameters. 
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